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ABSTRACT 

We present a comprehensive analysis of the relationship between star formation rate surface density, 
SsFR, and gas surface density, Sgas, at sub-kpc resolution in a sample of 18 nearby galaxies. We use 
high resolution Hi data from THINGS, CO data from HERACLES and BIMA SONG, 24 /xm data 
from the Spitzer Space Telescope, and UV data from GALEX. We target 7 spiral galaxies and 11 
late- type/dwarf galaxies and investigate how the star formation law differs between the H2-dominated 
centers of spiral galaxies, their Hi-dominated outskirts and the Hi-rich late- type/dwarf galaxies. We 
find that a Schmidt-type power law with index N — 1.0 ± 0.2 relates Ssfr and Sh2 across our 
sample of spiral galaxies, i.e., that H2 forms stars at a constant efficiency in spirals. The average 
molecular gas depletion time is 2 • 10^ years. The range of Sh2 over which we measure this relation 
is ~ 3 — 50 Mq pc^^, significantly lower than in starburst environments. We find the same results 
when performing a pixel-by-pixel analysis, averaging in radial bins, or when varying the star formation 
tracer used. We interpret the linear relation and constant depletion time as evidence that stars are 
forming in GMCs with approximately uniform properties and that Sh2 may be more a measure of 
the filling fraction of giant molecular clouds than changing conditions in the molecular gas. The 
relationship between total gas surface density (Sgas) and Ssfr varies dramatically among and within 
spiral galaxies. Most galaxies show little or no correlation between Shi and Ssfr. As a result, the 
star formation efficiency (SFE), SspR/Sgas, varies strongly across our sample and within individual 
galaxies. We show that this variation is systematic and consistent with the SFE being set by local 
environmental factors: in spirals the SFE is a clear function of radius, while the dwarf galaxies in our 
sample display SFEs similar to those found in the outer optical disks of the spirals. We attribute the 
similarity to common environments (low-density, low-metallicity, H i-dominated) and argue that shear 
(which is typically absent in dwarfs) cannot drive the SFE. In addition to a molecular Schmidt law, the 
other general feature of our sample is a sharp saturation of Hi surface densities at Ehi ~ 9 Mq pc~^ 
in both the spiral and dwarf galaxies. In the case of the spirals, we observe gas in excess of this limit 
to be molecular. 

Subject headings: radio lines: galaxies — radio lines: ISM — galaxies: ISM — galaxies: evolution 



1. INTRODUCTION 

A robust, quantitative measurement of the relationship 
between star formation rate (SFR) and gas density (SF 
law) is of major astrophysical importance in the context 
of galaxy evolution: it describes how efficiently galaxies 
turn their gas into stars and constrains theoretical mod- 
els of star formation; in addition, it serves as essential 
input to simulations and models of galaxy evolution 
(e.g ., Springcl & Hcrnquist 2003; Boissicr & Prantzoa 
19991: lTan et'al.1 Il999l : iKrumholz fc McKed I2005F 
Matteucci et al. r i2006D . Direct observations of this 



relationship at sub-kpc scales are still very rare because 
measuring the distributions of star formation, atomic 
gas, and molecular gas at high resolution and sensitivity 
is challenging and time-consuming. 

However, the past few years have seen an explosion 
in multiwavelength data for nearby galaxies. From 
the 'GALEX Nearby Galaxy Survey' (jGil de Paz et all 
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l2007t l and th e ^Spitzer Infrared Ne arby Galaxies Sur- 
vey' f SINGS. iKennicutt et aI1l2003f) the distribution of 
star formation is now known in a large suite of local 
galaxies. From the 'BIMA Survey o f Nearby Galax- 
ies' (BIMA SONG, iHelfer et all [200I ) and recent ob- 
servations with th e IRAM 30m telescope (HERACLES, 
iLerov et al.ll2008"H) the CO distributions of many nearby 
galaxies are also kn own. With 'The H I Nearby Galaxy 
burvey ' (THINGS, I Walter et al.l[200l . Hi maps that 
match or exceed the angular resolution and sensitivity 
of the ultraviolet, infrared, and CO data are now avail- 
able for 34 nearby galaxies. In this paper, we combine 
this suite of multiwavelength data to measure the sur- 
face densities of Hi, H2, and SFR (Shi, Sh2, and Ssfr) 
across the entire optical disks of 18 nearby galaxies at 
750 pc spatial resolution. Using these measurements, we 
examine the relationships among these three quantities 
across our sample. 

Following the pioneering work of lSchmidtJ (fl959h . it is 
common to relate gas density to the SFR density using a 
power law. He suggested the form psFR ~ (Pgas)", where 
PSFR and pgas denote the volume densities of the SFR and 
the gas. Studying the distribution of H i and stars per- 
pendicular to th e Galactic p l ane, h e der ived a pow er-law 
index of n « 2. ISanduleald (|1969( ) and iHartwicl (|1971l ) 
carried out the first measurements of the Schmidt law in 
other galaxies. They used bright stars in the Small Mag- 



2 



F. Bigiel et al. 



ellanic Cloud and Hii regions in M31, respectively, to 
trace star formation and focused their analyses on surface 
densities, Ehi and Esfr , which are directly observable 

(i.e., EsFR i^g.sf)- They found TVsMC = 1.84±0.14 
and Nm3i = 3.50 ± 0.12. Throughout this paper, we use 
the term Schmidt law or star formation law (SF law) to 
refer to a power law formalism relating star formation 
surface density to atomic, molecular, or total gas surface 
density. Note that for a constant scale height, the ex- 
ponents N and n are identical, i.e., it does not matter 
w hether one considers volume or surface densities. 

iMadore et al.1 ()1974D compared stars and Hll-regions 
to Ehi and derived a higher power law index for the o uter 
part than for the inner part of MSB. Newton. (jl980f ) re- 
peated this analysis wit h higher-resolution H i dat a and 
confirmed this behavior. iTosa &: Hamaiirnal ()1975l ) com- 
pared H ii-regions to Ehi in M31 and the LMC and find 
N 2 for both galaxies. iHamaiima fc Tosal (fl97l per- 
formed the same analysis in 7 nearby galaxies. They 
derived power law indices N — 1.5 — 2.9 and found radial 
variations in N for 2 of their sample galaxies (M31 and 
MlOl). 

iKennicuttI ()1989l . Il998a'. hereafter K98) studied the 
globally averaged relationship between SFR and gas in 
a sample of 61 nearby normal spiral and 36 infrared- 
selected starburst galaxies. K98 showed that a Schmidt 
law relates the disk-averaged total gas surface density, 
Egas = Ehi + Eh2, to the disk-averaged star formation 
surface density, Esfr, over many orders of magnitude. 
His subsample of normal spiral galaxies yields a power- 
law index N = 2.47 ±0.39; his composite sample, includ- 
ing starburst galaxies, yields N — 1.40 ± 0.15. Similar 
studies of the disk-averaged Schmidt Law used a range 
of SFR tracers - such as Ha, UV, radio cont inuum, and 
FIR emission - a nd found = 0.9- 1.7 (e.g.,|BiialeialJ 
[1989: Buat]fT99l iDeharveng eFallfiggl . 

Other authors studied a local Schmidt law using radial 
profiles (i.e., com paring azimuthally av eraged values) of 
EsFR and Egas. IWong fc Blitg (|2002[ ) fo und N = 1.2 



2.1 for 6 molecule-rich spiral galaxies. Boissie r et al 



( 2003 ) derived N « 2 for 16 galaxies and Misi riotis et al 



20061) obta i ned N = 2.18 ± 0.20 for the Milky Way. 
iHever et al.' ('2 00I f o und N = 3.3 ± 0.1 for M33. For 
2003) found TV = 1.4 ± 0.6 and 
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iKennicutt et all (|2007f L using 520 pc apertures centered 
on Ha and 24 /zm emission peaks, found N = 1.56±0.04. 

The large range of power-law indices in the literature, 
N Ki 1-3, suggests that either different SF laws exist in 
different galaxies or that N is very sensitive to system- 
atic differences in methodology (e.g., the choice of SFR 
tracers, the spatial resolution of the data, etc.). There 
are physically motivated, theoretical reasons to expect 
N in the range 0.75 - 2, and the pr ecise value of N may 
vary w ith the regime one considers. iKrumholz fc McKed 
(l2005h argue that a giant molecular cloud (GMC) will 
convert its gas into stars over a free-fall time, tg oc Pga°'^- 
If this is the case, then from the observati on that the sur- 
face density of GMCs is constant (e.g., I Solomon et al.l 
Il987t iBlitz et al.ll2007D . implying pgas oc M^^^, one 
expects N w 0.75 for the case where we measure the 
Schmidt law for individual GMCs. If we instead com- 
pare uniform populations of GMCs, where the different 
molecular gas surface densities reflect only different num- 



bers of clouds (not different physical properties) , then one 
expects N I for the molecular gas exponent (assum- 
ing the population-averaged timescale over which a GMC 
converts its gas into stars is constant). If gravitational 
instability in the neutral (H1+H2) ISM is the key pro- 
cess in star formation, the SFR may instead depend on 
the free-fall time of the total gas; in this case, and if the 
gas scale height is constant, tg oc Ega^-^ and one would 
expect N « 1.5 (;e.g.. lMador'ill977[ ). Finally, if one pos- 
tulates that star formation is a fundamentally coUisional 
process, because, e.g., it may depend on the formation of 
H2 by collisions between hydrogen atoms and dust grains 
or the collision of small clouds to form larger clouds, one 
would expect psFR oc Pgas, which will lead to w 2 
(again for a system with constant gas scale height). 

Another open question is whether the Schmidt law is 
fundamentally a molecular phenomenon or if a single, 
universal power law relates total gas and star formation. 
Because all stars are believed to form in molecular clouds, 
it would seem natural that H2 and SFR are more imme- 
diately related than Hi or total gas and SFR. There- 
fore it is somewhat surprising that observations remain 
contradictory on this point. One of the most surprising 
findings by K98 was that Esfr correlated better with 
Ehi than with Eh2 in normal disks; the strongest cor- 
relation, between Esfr and Egas, was o nly marginally 
strong er than the Esfr - Ehi correlation. IWong fc Blit j 
(|2002f ) found a much stronger relationship between Esfr 
and Eh2 than between Esfr and Ehi, even finding an 
anti-correlation between Esfr and Ehi at high SFRs, 
but they focused on molecule-rich spirals and therefore 
did not include large amounts of the Hi-dominated ISM. 
Even i n H2-rich spi rals, the questio n is still open. Re- 
cently , ISchuster et al. (2007) and Cr osthwaite fc Turneil 
()2007h studied the molecular gas-rich spirals NGC 5194 
and NGC 6946 and found that the total gas correlates 
better with the SFR than H2 alone. 

This paper attempts to answer the following questions: 
how do the SFR, H2, and Hi surface densities relate to 
one another in nearby galaxies on a pixel-by-pixel basis? 
Which of these relationships are common across our sam- 
ple, and which vary with environment? To address these 
questions, we measure Ehi, Eh2, and Esfr in a sample of 
18 galaxies: 7 spirals which have central areas (hereafter 
loosely referred to as centers) dominated by molecular 
gas ('spirals') and 11 Hi-dominated galaxies. We look at 
how these relationships differ between the H2-dominated 
centers of the spirals, their H i-dominated outskirts, and 
Hi-rich late type galaxies. By probing out to the optical 
radius r25, i.e., where the B-band magnitude drops below 
25magarcsec~^, and including late-type galaxies, we are 
able to strongly constrain the universality of the various 
Schmidt laws. 

Wc organize this analysis as follows: we describe the 
datasets used to measure Ehi , Eh2, and Esfr in §[21 
We explain how we convert these data to physical units 
and how we generate a set of independent measurements 
over the disk of each galaxy. In §[3] we show the observed 
relationships between Ehi, Eh2, and Esfr for individual 
galaxies. We carry out power law fits to each distribu- 
tion and assess whether there is a common Schmidt law 
shared by all of our galaxies. We show results for spi- 
rals and Hi-dominated late- type galaxies separately to 
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highlight the range of environments in our sample. In 
§|4]we combine our data for the spirals and for the Hl- 
dominated galaxies, to be able to draw general conclu- 
sions about Hi, H2, and the SFR in these two samples. 
In §[5] we explore the effect of environment by compar- 
ing the Schmidt law in the inner, H2-dominated parts 
of spirals to their H i-dominated outer disks and to H i- 
dominated dwarf and late- type galaxies. We show that 
the relationship between total gas and SFR is a clear 
function of environment and fit the radial dependence of 
the molecular-to-atomic gas ratio in spirals. In §|6] we 
summarize our results and give our conclusions. Here we 
also compare our results to those from previous studies. 

2. SAMPLE, DATA, UNITS, AND MEASUREMENTS 

This work is made possible by the existence of the 
state-of-the-art multiwavelength data sets mentioned in 
the introduction. We use high resolution Hi data from 
'The Hi Nearby Galaxy Survey' fTHINGS- IWalter et all 
I2OO8I ) and CO maps from a large ongoing project us- 
ing HERA on the IRAM 30m telescope (HERACLES, 
iLerov et all l2008bD and from the BIMA SONG survey 
(BIMA SONG. Heifer et al.l[200l . We trace star forma- 
tion combining mid-IR maps fro m the 'Spitzer Infrared 
Nearby Galaxies Survey' f SINGS. [Kennicutt et all2003f) 
and FUV n iaps from the 'GALEX Nearby Galaxy Sur- 
vey' (NGS, iGil de Paz et allfeOOTl ). In this section we 
describe our sample and these datasets, show maps and 
radial distributions, and explain how we generate mea- 
surements of Shi, 2]h2, and Ssfr- 

2.1. Sample 

Our sample of spiral galaxies consists of the cross- 
section of THINGS, SINGS, the GALEX NGS and the 
HERACLES or BIMA-SONG CO surveys. Galaxies that 
are known to be H i-dominated are not required to be 
part of HERACLES or BIMA-SONG. This ensures that 
we measure at least the dominant component of the neu- 
tral ISM and the SFR along most lines of sight. Table 
[1] lists our sample along with their adopted pro perties: 
distance, inclination, position angle, radius f,Walte r et al.l 
I2OO8I except that we adopt i = 20° in M51) and morphol- 
ogy (from LEDA. IPrugniel fc Heraudeaull 19981) . We sep- 
arate the galaxies into two groups: 7 large spiral galaxies 
that have an H2-dominated ISM in their centers ('spi- 
rals') and 11 late-type galaxies that are Hi-dominated 
throughout. 

We do not study galaxies with inclinations > 70° that 
would otherwise qualify for the sample. High inclinations 
yield very few sampling points and result in a deprojected 
beam elongated parallel to the minor axis in the plane 
of the galaxy, making it difficult to carry out fits to the 
data and interpret the results. We also do not include any 
galaxies more distant than 12 Mpc. This is the maximum 
distance at which the typical angular resolution of our 
data corresponds to our common spatial resolution of 
750 pc. 

We work without CO maps for 4 late-type spirals, 
NGC 925, NGC 2403, NGC 2976, NGC 4214 and 7 dwarf 
irregular galaxies. For most of these galaxies, the CO 
content is either measured or constrained by a signif- 
icant upper limit: major axis profiles by lYoung et al.l 
(fl995[ ) for NGC 925, NGC 2403, and NGC 2976 (which 
is also part of BIMA SONG); Kitt Peak 12m mapping of 



TABLE 1 

Sample Properties-'- 



Galaxy D i PA 

[Mpc] [dcg] [dcg] 



r25 Hubble 
[arcmin] type 



H I-dominated Galaxies 



DD0154 


4.3 


66 


230 


0.97 


Irr 


DD0 53 


3.6 


31 


132 


0.39 


Irr 


Hoi 


3.8 


12 


50 


1.66 


Irr 


HoH 


3.4 


41 


177 


3.30 


Irr 


IC 2574 


4.0 


53 


56 


6.44 


SABm 


M 81 DwA 


3.6 


23 


49 


0.63 


Irr 


M 81 DwB 


5.3 


44 


321 


0.56 


Irr 


NGC 925 


9.2 


66 


287 


5.36 


Scd 


NGC 2403 


3.2 


63 


124 


7.92 


SABc 


NGC 2976 


3.6 


65 


335 


3.62 


Sc 


NGC 4214 


2.9 


44 


65 


3.38 


Irr 


Spirals With H2-Dominated Centers 


NGC 628 


7.3 


7 


20 


4.89 


Sc 


NGC 3184 


11.1 


16 


179 


3.71 


SABc 


NGC 3521 


10.7 


73 


340 


4.16 


SABb 


NGC 4736 


4.7 


41 


296 


3.88 


Sab 


NGC 5055 


10.1 


59 


102 


5.87 


Sbc 


NGC 5194 


8.0 


20 


172 


3.88 


Sbc 


NGC 6946 


5.9 


33 


243 


5.74 


SABc 



^ See Walte r et al.l l|2008l ^ for further information on in- 
dividual galaxies and for references to the values quoted 
in this table. 



the inner disk of NGC 2403 (jThornlevfc Wilsonl[l99| ): 
OVRO mapping of NGC 4214 bv IWalter et all (|2f)nih : 
single dish point ings toward NGC 4214, Holmberg I 
and DDO 154 bv iTavlor et~an (|1998D : centr al pointings 
forDDO 53, IC 2574, and Holmberg II by ILerov et al.l 
(|2005h . In each case, the ISM is well-established to be 
H I-dominated. M81 DwA and M81 DwB, two extremely 
low- mass dwarf irregular companions to M81, both lack 
CO measurements, but should also be Hi-dominated. 
Of these galaxies, 8 have absolute _B-band magnitudes 
fainter than —18 mag and maximum rotation velocities 
< 100 km s^^ (i.e., they have the mass of the LMC or 

less). NGC 925 and NGC 2403 have Mb 20 mag 

and maximum rotational velocities ~ 100 km s^^; they 
are probably intermediate in mass between the LMC and 
M 33. 

2.2. Data 
2.2.1. THINGS Ri Maps 

To estimate the surface density of neutral atomic hy- 
drogen, Shi, we use VLA maps of the 21 cm line obtained 
as part of 'The H i Nearby Galaxy Survey' (THINGS, 
IWalter et ani2008D . THINGS consists of high-resolution, 
high-sensitivity H i data for 34 nearby galaxies obtained 
with the NRAO*" VLA. The target galaxies have dis- 
tances of 2 < D < 15 Mpc and cover a wide range in 
star formation rates, H i masses, luminosities and Hub- 
ble Types. 

The FWHM of the primary beam (field-of-view) is 
32'. We use 'robust' weighted maps, which have a typ- 
ical beam size of '^6". The typical la RMS noise at 
our working resolution of 750 pc (see ^ \2.'S^ is cr(SHi) « 
0.5 M0 pc~^. Because our analysis is restricted to the 
regime within the optical radius r25 (see ii l2.4l and com- 
pare the radial profiles in ii l2.5p . the measured Hi sur- 



° The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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face densities are safely above this sensitivity limit for 
all galaxies. 

The integrated Hi intensity map for NGC 6946 at 
our working resolution of 750 pc is shown as an exam- 
ple in the left panel of Figure [TJ Further details regard- 
ing the data produ cts and data reduction are given in 
IWalter et alJ (|200l . 

2.2.2. HERACLES CO Maps 

We derive Sh2 distributions from two sources: an 
ongoing nearby galaxy survey using the HERA focal 
plane array on the IRAM 30m telescope (HERACLES, 
iLeroY et all l2008bf) . and the 'BIMA Surv ey of Nearby 
Galaxies' (BIMA SONG. lHelfer et al.l[2003h . The profiles 
from the two sets of CO maps agree with one another and 
also with the major axi s pointings from the FCRAO sur- 
vey (|Young et al.lll995D . We use the HERACLES maps 
where available because they have good extent and sen- 
sitivity, often measuring CO well into the H i-dominated 
regime where (on average) Efji > Sh2 7 i-e., beyond the 
H2-to-Hi transition radius. Extending beyond the tran- 
sition radius is critical to be able to differentiate between 
a molecular and total gas Schmidt law. 

We use HERACLES maps for 6 spiral galaxies: 
NGC 628, NGC 3184, NGC 3521, NGC 4736, NGC 5055 
and NGC 6946. These were obtained at the IRAM 30m 
telescope during Januar y and Octob e r 2007 and further 
details are described in iLeroy et al.l ()2008bl ). We show 
the CO intensity map for NGC 6946 in the middle panel 
of Figure [TJ HERA maps the CO J = 2 ^ 1 transition 
with an angul ar resolution o f 11" . For technical details 
on HERA see ISchuster et al.) ()2004l ). Our maps are sen- 
sitive to surface densities Eh2 ~ 3 M© pc~^. 

BIMA SONG produced maps of CO J = 1 ^ emis- 
sion with good resolution, ~ 7", but limited field of view. 
We use these maps only for NGC 5194. Also, we use 
these data to compare CO to the SFR at high resolu- 
tion in the central parts of all 7 centrally H2-dominated 
galaxies. The BIMA SONG map for NGC 5194 used 
here includes zero-spacing data from the Kitt Peak 12m 
telescope, making it sensitive to extended structure. The 
BIMA SONG maps we use are sensitive to surface den- 
sities above Eh2 ~ 10 M0 pc~^. 

2.2.3. GALEX FUV Data 

We use FU V data from the 'GALE X Nearby Galaxy 
Surve y' (NG S, iGil de Paz et al.ll2007f) to estimate Ssfr 
(see ii l2.3.2p . The GALEX observatory provides simulta- 
neous imaging in a far UV (FUV) and a near UV (NUV) 
broadband filter. The FUV band covers the wavelength 
range 1350- 1750 A, the NUV covers 1750- 2800 A. The 
angular resolutions (FWHM) are 4.0" and 5.6" for the 
FUV and the NUV respectively, and the field-of-view of 
the instru ment is 1.25°. For tech nical details on the de- 
tector see iMorrissev et all (|2005D . 

We remove a single set of foreground stars from the 
GALEX and Spitzer 24/im maps (see ^ I2.2.4[) . Fore- 
ground stars are prominent in the NUV images but much 
less so in the FUV images, and we find that we can iden- 
tify them easily via their UV colors. We therefore use 
the NUV maps to identify foreground star candidates 
and apply an NUV/FUV ratio cut of NUV/ FUV > 10. 
We apply this to all emission above 5 a in the NUV maps, 
below this we cannot distinguish stars from noise. In a 



few cases, a cut-off higher than 10 was necessary to re- 
move particularly bright foreground stars. We estimate 
the contribution from foreground stars below 5 cr in our 
FUV images beyond 2 r25 (where there is only marginal 
emission from our target galaxies) and find that their 
contribution is negligible. 

We estimate and remove a small background from the 
FUV maps. We measure this away from the galaxy by 
discarding emission with intensities > 3a above the me- 
dian value of the image. The background adopted is the 
median of this residual map and is subsequently sub- 
tracted from the original map. The backgrounds in the 
GALEX maps are very well behaved and this simple sub- 
traction yields good results in most cases. For three 
galaxies - NGC 4214, NGC 5194 (M51), and NGC 6946 
- we blank portions of the map that show the edge of 
the GALEX ficld-of-view, obvious artifacts, bright stars 
not entirely removed by the color-cut as well as M51b. 
These are usually located well away from the galaxy and 
have minimal impact on this study (which is restricted 
to the optical disks of the galaxies). 

We correct the FUV maps for the effects of Galactic 
extinctio n. We estimate E(B — V) from the extinction 
maps of ISchlegel. Finkbeiner fc Davis! (jl998D and con- 
vert to FUV extinction using Apuv = 8.24 x E{B - V) 
(Wyder et al. 2007). 

2.2.4. Spitzer Space Telescope 24 fJ^m Data 

FUV data are heavily affected by internal extinc- 
tion. To estimate the amount of star formation obscured 
by dust (see S I2.3.2p . we use 24 ^m data obtained by 
the ^Spitzer Infrared Nearby Galaxies Survey' (SINGS, 
iKennicutt eraI]|2003D . The map of NGC 4214 comes 
from the public archive. These data are scan maps taken 
with the MIPS instrument o n board the Spitzer Space 
Telescope (|Rieke et al.ll2004D . Full details of the ob- 
serving strategy and t arget sensitivity are described by 
IKennicutt et al.l ()2003[ ). The MIPS data were processed 
using the MIPS Ins trument Team Data Analysis Tool 
(|Gordon et al.l[200l . 

The SINGS observations were designed to detect emis- 
sion out to the optical radius and in most of our targets 
this goal was achieved. The FWHM of the MIPS PSF 
at 24 /im is 6", though the beam is substantially non- 
Gaussian outside the central peak. This is only a minor 
concern for this study because to achieve a common spa- 
tial resolution for all data, we work with a typical angular 
resolution of ^ 20". 

We perform the following additional processing on the 
SINGS maps: we first blank the edges of the scan maps 
parallel to the direction of the scan; in these regions the 
noise increases and artifacts are more common. We also 
blank the stars identified from the GALEX images. We 
then subtract a background, which is usually quite small, 
using the same procedure as for the GALEX images. In a 
few maps - Ho I, NGC 3521 and NGC 6946 - we identify 
regions with unreliable background and blank these by 
hand. These regions are far away from the galaxy though 
and so have minimal effect on this analysis. 

2.3. Alignment, Units, and Convolution 

From the data described above, we construct maps of 
star formation rate surface density, Ssfr, atomic hy- 



The Star Formation Law in Nearby Galaxies on Sub-Kpc Scales 



5 



drogen surface density, Shi, molecular hydrogen surface 
density, SH2, and total gas surface density, Sgas- 

We place all of these maps on the THINGS astromet- 
ric grid (pixel scale: 1.5") at a common spatial reso- 
lution of 750 pc. We carry out most of our work at 
750 pc resolution but also create versions of each map at 
a range of spatial resolutions, from the native resolution 
to 10 kpc (in steps of 50 pc below 1 kpc, and 500 pc up 
to 10 kpc). When degrading the resolution of a map, we 
convolve with a circular Gaussian beam (on the sky), i.e., 
we do not account for the inclination of the galaxy. For 
a given spatial resolution, this exercise may be thought 
of as moving each galaxy at a larger distance. 

2.3.1. Gas Surface Density Maps 

The total gas surface density, Sgas, is the sum of the 
atomic gas surface density, Shi, and the molecular gas 
surface density, Sh2 (we do not consider ionized gas, 
as this makes up only a small fraction of the total gas 
content). Where Sh2 is below the sensitivity limit (see 
§ I2.2.2|) . we take Sgas = Shi, which is formally a lower 
limit. These quantities all have units of Mq pc~^ and 
are hydrogen surface densities, i.e., they do not include 
any contribution from helium. To scale our quoted sur- 
face densities to account for helium, one should multiply 
them by a factor of ^ 1.36. All surface and column den- 
sities quoted in this paper have been corrected for the 
inclinations given in Table [TJ 

We assume a ratio I {2 1)/^(1 ^ 0) = 0.8, a typ- 
ical num ber derived fr om c o mparing our maps t o the 
maps of Eelfer et all (|2003f l. lYoung et all (| 19951 ) and 
iKuno et"al] (|2007f ). A det ailed compar i son bet ween these 
datasets can be found in iLerov et al.l (|2008bD . To con- 
vert from Sco to Sh2, we adopt a CO-to-H2 conversion 
factor of 2.0 x lO^^ cm^^ (j^ j^j^ s-^)-'^. This value is 
appropriate for th e Milky Way according t o 7-ray and 
FIR studies (e.g., IStrong fc Mattoxl 119961 : iDame et all 
|2001| ). but slightly lower than the c onversion factor sug- 
gested by virial m ass methods (e.g.. lSolomon et al.lll987l : 
iBlitz et al]|2007[ ). We do not account for any changes 
in the conversion factor with metallicity or other en- 
vironmental factors. For comparison, we note that we 
use the same conversion factor adopted bv lWong fc Blit j 
(|2b02i) but that K98 adopted a higher value of 2.8 x 
10^0 cm-2 (K km s-^)-^. 

2.3.2. Star Formation Rate Surface Density Maps 

We measure star formation rate surface densities, 
SsFR, in units of M0 yr^-'^ kpc^^, by combining FUV 
and 24/xm maps. The GALEX FUV images trace pre- 
dominantly O and early B stars and therefore offer a 
picture of recent, unobscured star formation. FUV emis- 
sion can be heavily affected by dust, however, making it 
difficult to derive unbiased measurements of Ssfr from 
the FUV alone. Fortunately, the dust that obscures the 
FUV emission from young stars is heated and reradiates 
in the mid-infrared. The 24 fim maps thus allow us to 
estimate the amount of ongoing dust-obscured star for- 
mation. 

This approach was proposed and validated for indi- 
vidual star formin g regions and la rge p ortions of galac- 
tic dis ks by iCalzetti et al.l (|2007( ) and iKennicutt et al.l 
(|2007f ). They showed that 24 /xm emission could be used 



to accurately estimate the amount of extinction that af- 
fects Ha emission. We adopt an analogous approach and 
combine FUV and 24 /im using the following formula: 

Ssfr [M© yr'^ kpc"^] =3.2 x 10^^/24 [MJy ster"^] (1) 

+ 8.1 X IQ-^Ipuv [MJy ster"^] , 

where Ssfr denotes the star formation rate surface den- 
sity, and I24 and Ifjiv are the 24 /im and FUV intensities 
respectively. iLerov et al.l (l2008a) motivate the choice of 
coefficients in Equation [T] Ssfr derived from Equation 
[1] agrees with th e values derived usin g the Ha -I- 24 /im 
calibration from ICalzetti et al.l (|2007f ) at 750 pc resolu- 
tion; when /24 = 0, Equati on [1] reduces to the FUV-SFR 
calibration bv lSalim et al.l ([2007). SFRs integrated over 
entire galaxies as well as values derived from azimuthally 
averaged radial profiles agree well with those based on 
extinction-corrected Ha emission alone. 

As an example we show the Ssfr map for NGC 6946, 
which was derived using Equation [1] in the right panel 
of Figure [TJ The RMS noise varies from map to map, 
but a typical limit for the sensitivity of the SFR maps is 
« 1 X lO"'' M0 yr-i kpc" ^ For our calib r ation of Ssfr 
we adopt the IMF from ICalzetti et"all (120071), which 
is th e default IMF in STARBURST99 ([Leitherer et all 
fT999h . i.e., a Kroupa-type two-component IMF t hat ex- 
tends to 120 M0. To conver t to the truncated iSalpeted 
(1955 ^ IMF adopted by , e.g.. lKennicufrg (|1989lll998af) or 
iKennicutt et al.l (|2007| ). one should multiply our Ssfr 
by a factor of 1.59. 

Because this particular combination of FUV and mid- 
IR maps is new, we compare our results to those ob- 
tained using maps of Ha and Ha-|-24/im emission wher- 
ever these are available. For this purpose, we use Ha 
maps from the SINGS Data Release 4. We fit and re- 
move backgrounds from these maps, masking out any 
region with a particularly problematic ba ckground. We 
correc t for [Nil] in the bandpass following ICalzetti et al.l 
(l2007f) and for Galactic extinc t ion us ing the dust maps of 
ISchlegel. Finkbeiner fc David (|1998l ). We use Ha maps 
for all of our spiral galaxies except NGC 4736. 

In addition to sampling our data on a pixel-by-pixel 
basis at 750 pc resolution, we also show the results of 
our analysis carried out in azimuthally-averaged radial 
profiles and explore the efi^ects of degrading to lower res- 
olution (so that each 'pixel' contains many square kpc). 
Further, we compare our results to measurements from 
other studies using a variety of SF tracers (see ii l6.2p . 
These tests are shown in the appropriate sections below 
and lead us to conclude that our SFR maps are accu- 
rately tracing the amount and distribution of recent star 
formation. 

2.4. Individual Data Points 

With maps of Shi, Sh2, and Ssfr in hand, the next 
step is to generate individual data points. We com- 
pare these quantities over the entire optical disk of each 
galaxy, out to r25 — ^25/2. This is a departure from pre- 
vious studies, which used integrated measurements over 
entire galaxy disks (e.g., Kcnnicutt 1989, 1998a,), radial 
profiles (e.g.. lWong fc Blit^l2002 ). or a pertures centered 
on H a and 24 /im emission peaks (e.g., IKennicutt et al.l 
I200I . 
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Fig. 1. — Maps of integrated Hi (left), CO (middle) and SFR (right) surface densities at our working resolution of 750 pc for the spiral 
galaxy NGC 6946. The ellipses indicate the optical radius (r25) in the plane of the galaxy, within which we compare Sjji, Sh2i and Sgpjj^. 
Almost all star formati on occurs within this radius, although the Hi often extends much beyond r25. The maps for all other galaxies in 
our study are shown in lLerov et al.l l l2008a ). 



We choose r25 as an outer limit. Most star formation 
takes place within t his ra dius (e.g., see the p rofiles of 
iMartin fc KennicuttI [200l IWong fc Blitzl l2002t ) and our 
FUV, 24/im, and H i maps detect emission at more than 
3(7 over most of this area. The ellipses in Figure [T] show 
the projected optical radius, i.e., the extent of the re- 
gion that we study, on the THINGS H i, HERACLES CO 
and SsFR maps for NGC 6946. We refer the reader to 
iLerov etaTI (|2008al) for a 'galaxy atlas' showing all Shi, 
Eh2, and Ssfr maps for the galaxies in our sample. 

We draw independent data points from our maps so 
that each data point corresponds to non-overlapping res- 
olution elements and the data together cover the opti- 
cal disk. This may be thought of as either reducing the 
oversampled map (rebinning the map) so that one pixel 
corresponds to a resolution element or as covering the 
disk with non-overlapping apertures equal in size to a 
resolution element. For each independent data point we 
measure Shi, Sh2, Sgas = Shi + Sh2, and Ssfr- 

When convolving the maps with progressively larger 
beams (see ^ \2.3\\ , flux from outside our radius cut at 7-25 
would be convolved into the optical disk. To avoid this, 
the radius cut we use is decreased by half a beam width 
and is actually rmax = ^25 - 0.5 6'boam, where Obeam is 
the FWHM of the beam in the same units as the radius. 
Note that ^beam is subtracted in the plane of the sky, not 
the plane of the galaxy because the convolution takes 
place in the plane of the sky. 

2.5. Radial Profiles 

We also extract from our maps radial profiles of Shi, 
Sh2, and Ssfr and examine them along with the pixel- 
by-pixel data. We show the profiles for the spirals and 
8 of our H i-dominated galaxies in Figures [2] and [31 We 
note that the radial profiles as well as the Hi, CO and 
SFR m aps for all of our galaxies are shown in iLerov et al.l 
(|2008bf) . The plots in Figures [D and [3] show Ssfr, Shi 
and Sh2 as a function of galactocentric radius normalized 
by ^25. The top axis gives the radius in kpc. Each point 
represents the average value in an individual 10"-wide 
tilted ring within 60°of the major axis using the struc- 
ture parameters given in Table [T] For the Hi-dominated 



galaxies, we show only Ssfr and Shi. 

3. THE STAR FORMATION LAW IN INDIVIDUAL 
GALAXIES 

Figure [5] shows the relationship between gas and star 
formation surface densities in individual spiral galaxies. 
Color contours show the results from pixel-by-pixel sam- 
pling the optical disks of the 7 spiral galaxies in our 
sample. Each row shows results for a different galaxy. 
The columns left to right show Ssfr on the y-axis as a 
function of Shi (left), Sh2 (middle), and Sgas (right) on 
the X-axis. Shading shows the density of data points in 
cells of size 0.05 dex wide (in both axes). Green, orange, 
red, and magenta cells show contours of 1, 2, 5, and 10 
sampling points per cell. Individual cells that are not 
connected and contain only one data point are marked 
in grey. All data in Figure 3] are at a common resolution 
of 750 pc. 

We plot points from the radial profiles (Figures [2] and 
[3]) on the same plot as black crosses. Generally speak- 
ing, the radial profile data lie near the middle of the 
distribution of the pixel-by-pixel data, as expected. In 
many cases, features at high Sgas and high Ssfr are 
clearer in the radial profile data because small central 
rings with high surface densities have small areas and 
thus contribute only a few individual points to the pixel- 
by-pixel data. 

A common way to parametrize the relationship be- 
tween star formation and gas is via the gas depiction 
time or its inverse, the star formation efficiency SFE = 
SsFR/5^gas- The gas depletion time is the time needed for 
the present rate star formation to consume the existing 
gas reservoir. In Figure 31 we plot levels of constant de- 
pletion time/SFE as diagonal dotted lines. From bottom 
to top, these correspond to gas depletion times of 10^", 
10^ and 10^ years or equivalently depleting 1%, 10% and 
100% of the gas reservoir per 10® years; these values in- 
clude helium and so are true depletion times. 

3.1. Fits to the Ssfr Versus Sgas Distributions 

Solid black lines in Figure 2] show the results of fitting 
a power-law of the form 
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Fig. 2. — Azimuthally averaged radial profiles of Shi, ^H2, and SgpR for spiral galaxies with H2-dominated centers. The y-axis shows 
Shi and Sh2 in units of Mq pc~^ as well as SgpR in units of 10~^ Mq yr~ kpc~^ (the scaling is chosen to bring the profiles onto the 
same plot). The s-axis shows galactocentric radius normalized by ^25 (bottom) and in kpc (top). Profiles of Ssfr and Sh2 are strongly 
covariant. Shi varies weaJjly over the optical disk with the main features often being a central depression and a universal upper limit of 
Shi » 9 Mq pc'^. 
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Fig. 3. — Azimuthally averaged radial profiles of Shi and Ssfr iov 8 of our 11 H l-dominated galaxies, all late-type spirals or dwarf 
irregulars. The j;-axis shows Ejji in units of Mq pc^'^ or Sgpji in units of 10^'* Mq yr^-*^ kpc^'^ (the scaling is chosen to bring the 
profiles onto the same plot). The x-axis shows the galactocentric radius normalized by r25 (bottom) and in kpc (top). Shi shows the same 
maximum value seen in spiral galaxies and similarly flat profiles. 
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Fig. 4. — SgpR as a function of Shi (left), Sh2 (middle), and Sgas = Sjji + Sh2 (right) in our spiral galaxies at 750 pc resolution. 
Each row shows results for one galaxy. Green, orange, red, and magenta cells show contours of 1, 2, 5, and 10 independent data points 
per 0.05 dex-wide cell (for H2 in NGC 4736 we use a scatter plot due to the low number of sampling points.). Crosses show average 
measurements over tilted rings from the radial profiles. Diagonal dotted lines show lines of constant SFE, indicating the level of Sgpj^ 
needed to consume 1%, 10% and 100% of the gas reservoir (including helium) in 10* years. Thus, the lines also correspond to constant gas 
depletion times of, from top to bottom, 10*, 10^ , and 10^'^ yr. Dashed vertical lines in the Hi (left) and total gas (right) plots show the 
surface density where the Hi saturates (see S lS.Sj l. Dotted vertical lines in the middle plots show the typical sensitivity for our CO data. 
We show OLS bisector fits to the H2 and total gas data with a solid line and quote the results. 
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Fig. [H — continued. 
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using the ordinary least-squares (OLS) bisector. The 
free parameters are the power law index, TV, and a, which 
is SsFR at the fiducial gas surface density of 10 Mq pc^^. 
Figure H] and Table [2] give the best fit values of N ^ 
A — logjQ a, and the RMS scatter in log Ssfr about the 
fit. Because there is not a clear independe nt variable, 
we ca rry out the fit using the OLS bisector (jlsobe et al.l 
Il990l ). giving equal weight to each point. We treat the 
problem in log space, where fitting a line yields N as 
the slope and A as the intercept. By centering the fit at 
10 M0 pc~^, a surface density near the middle of the dis- 
tributions, we minimize the covariance between N and 
a. We note that readers interested in comparison with 
previous work should take note of this difference; most 
fits in the literature quote 'A' at 1 M© pc~^ rather than 
our fiducial 10 Mq pc^^. 

Wc fit SsFR vs. Sh2 (middle panel) in the regime 
Eh2 > 3 Mq pc-2 for the HERACLES data and Sh2 > 
10 MoPC-2 for the BIMA SONG data (see ^ [2X2]) . 
These sensitivities are shown as a dotted vertical line 
in the middle panels of Figure ID We fit Esfr vs. Egas 
to all data points (right panel) . The formal errors on the 
fits are small and the fits are robust to the removal of 
individual data points. We test the latter via bootstrap- 
ping (i.e., repeatedly drawing a new, equal-sized, random 
subsample from our data allowing repetition) and find 
the resulting uncertainties to be typically ~ 0.01 in N 
and ~ 0.05 in A. Methodology, e.g., the decision to fit 
X vs. Y, Y vs. X, or use of the OLS bisector, drives 
the resulting fits as much as any other statistical factor. 
Because both variables are independent, we use the OLS 
bisector, but estimate the uncertainty by carrying out or- 
dinary least-squares (OLS) fits of Ssfr vs. Sh and Sh 
vs. SsFR. The differences in TV and A that we obtain are 
taken as the uncertainties on the fit. We find that this 
brackets the range of reasonable 'by eye' fits well. We 
note that these uncertainties represent only the uncer- 
tainty in fitting the distribution of data points. They do 



not refiect systematics such as uncertainty or variations 
in the CO-to-H2 conversion factor, the IMF, etc. 

3.2. The Molecular Gas Schmidt Law 

Fitting SsFR to Eh2 alone yields power law indices 
near unity, N = 0.96 ±0.07 and coefficients A = -2.06 ± 
0.17 with a typical scatter of ~ 0.2 dex. That is our 
7 spirals display power law indices consistent with an 
= 1 molecular Schmidt law and only mild variations 
in the normalization. Another statement of this is that 
the molecular gas in our sample shows a nearly constant 
ratio of Esfr to Eh2 that corresponds to star formation 
consuming the H2 gas reservoir in ^ 2 x 10^ yrs (i.e., a 
constant gas depletion time; also see iS I4.3|) . 



3.3. The Total Gas Schmidt Law 

The relationship between Egas and Esfr shows a much 
larger range of behavior as in the H2 case. The power law 
indices, A^, for these fits range from 1.11 to 2.74 with the 
mean 1.85 ± 0.70 and a scatter of ^ 0.3 dex. This range 
is similar to that found by studies of individual galaxies 
in the literature (1 < A^ ^ 3, see §[1]). These steeper 
fits, as compared to the H2 case, are mainly caused by a 
drop in EgFR over a relatively narrow range of gas surface 
densities just below Egas ~ 9 Mq pc~^ (see §[33] where 
we assess this further ). 

There is less variation in the coefficient. A, than in the 
power law index. A varies from —2.75 to —2.00 with a 
mean of —2.39 and an RMS scatter of 0.28. This mean 
value is lower than the mean value for the molecular gas 
A = —2.06 (see i) l3.2p . The ratio of mean coefficients 
a for the total gas and the H2 is ~ 0.3 dex; that is, a 
particular surface density of molecular gas will on average 
form stars at about twice the rate of the same surface 
density of total gas. 

It is further evident that there is no universal behavior 
in the right hand column of Figure H) That is, the dis- 
tribution of points in Egas-EgFR parameter space varies 
from galaxy to galaxy. In some cases, e.g., NGC 5194, 
a single power law appears to relate the two. In other 
cases, e.g., NGC 3184, Esfr and Egas are essentially un- 
correlated where Ehi > Eh2. In these cases, there is no 
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TABLE 2 

Fitted Power-Law Parameters at 750 pc Resolution 



Ha Hi + Ha 

Galaxy Coefficient (A) Index (A'') Scatter Coefficient (A) Index (N) Scatter 

NGC628 -1.99 0.98 0.16 -2.35 2.74 0.39 

NGC3184 -2.16 1.12 0.18 -2.45 2.50 0.31 

NGC3521 -2.19 0.95 0.10 -2.75 2.12 0.19 

NGC4736 -1.79 0.95 0.15 -2.00 1.44 0.19 

NGC5055 -2.22 0.92 0.10 -2.63 1.58 0.22 

NGC5194 -2.09 0.84 0.19 -2.22 1.11 0.25 

NGC6946 -1.94 0.92 0.23 -2.33 1.46 0.29 



Average -2.06 ±0.17 0.96 ili 0.07 ■ ■■ -2.39 ili 0.28 1.85 ±0.70 



clear one-to-one relationship between total gas and SFR 
across the whole disk. 

What causes some galaxies to display a power law 
spanning from the H i-dominated to H2-dominated ISM 
while others do not? The data-driven answer is that 
galaxies which show a single power law have compara- 
tively low Hi surface densities within their optical disks. 
Galaxies with uncorrelated Sgas-SsFR distributions (e.g., 
NGC 3184) by contrast do not have many lines of sight 
with H I at low surface densities within their optical disks. 

One possible explanation is that galaxies with well- 
defined total gas Schmidt laws (and low H I surface den- 
sities) may have lost diffuse Hi unassociated with star 
formation in interactions. Many show signs of tidal dis- 
ruption or ongoing interactions. Assessing the underly- 
ing reason for the range of distributions in the right hand 
part of Figure H] is beyond the scope of this paper, the 
key observation here is that there is a range of distribu- 
tions and that many galaxies are not well described by a 
single power law relating Sgas to Ssfr- 

3.4. Star Formation Efficiencies 

The SFE provides another way to express the results 
discussed above. Most galaxies show a fixed SFE relat- 
ing their Ssfr and Eh2- Some galaxies (e.g., NGC 5194) 
also show a constant SFE in their total gas spanning from 
high to low surface densities. The case of NGC 5194 is 
particularly striking; this galaxy displays a nearly con- 
stant SFE spanning surface densities from 1 Mq pc~^ to 
100 Mq pc~^. Over two orders of magnitude in gas sur- 
face density, the power law index remains near unity and 
Ssfr shows only a factor of '--^ 2 scatter about a constant 
SFE. 

On the other hand, several galaxies, e.g., NGC 628, 
NGC 3184 or NGC 3521, show steep distributions in the 
right hand column of Figure [4] and high power law in- 
dices. This may be phrased as large internal variations 
in their SFE. The variations in SFE are as striking in 
their own way as the correlation in NGC 5194: these 
galaxies span nearly an order of magnitude in SFE at an 
almost constant Egas ~ 5 Mq pc~^. This is clear ev- 
idence that the total gas surface density cannot be the 
critical quantity setting the SFR over the H i-dominated 
parts of these galaxies. 

There is also variation in the SFE among galaxies. This 
can be seen from the range of coefficients to our power 
law fits. At a particular Egas, the average Esfr shows an 
RMS scatter of ~ 0.3 dex. Galaxy-to-galaxy variations 
thus account for a factor of ~ 2 scatter in the SFE in our 
sample. 



3.5. H I Saturation at High Column Densities 

Figure [4] and the radial profiles in Figure [2] also illumi- 
nate the relationship between Ehi and Eh2- Both show 
a striking absence of high surface density H i; this is seen 
as a sharp right-hand edge to the distributions shown in 
the left-hand column of Figure S] and the failure of our 
radial profiles (spirals or H i-dominated galaxies) to cross 
Ehi ~ 9 M0 pc~^. The only gas in excess of this limit- 
ing surface density a. ppears to be in the molecular phase. 
IWong fc Blitd (|2002D showed a similar 'saturation' effect 
in azimuthally averaged profiles for their molecular gas- 
rich spirals, as did, e.g., Martin fc Kennic utt (2001) and 
iMorris fc Lol (|1978[ ). The data plotted in Figure H show 
that this effect is present at 750 pc resolution and that 
it is remarkably universal. Esfr and Eh2 show no com- 
parable limiting values. 

The vertical dashed line in the left and the right 
columns of Figure H shows EH/,saturation ~ 9 Mq pc"^. 
In Section § 14.11 we will see that 95 % of the Hhi values 
for the combined distribution of all 7 spiral galaxies in 
our sample are below EH/,saturation- 

A second effect is best seen in the radial profile points 
in the left column of Figure SI at high SFRs there is 
often an anii-correlation between Shi and Ssfr- This 
occurs in the central H l-holes of spirals where the gas is 
overwhelmingly molecular and the star formation rate is 
very high. All 7 of our spiral galaxies show some degree 
of this effect, i.e., at least a mild central depression in 
Hi. 

3.6. Ssfr Versus Sgas in H i-dominated Galaxies 

We have already seen that there are variations in the 
relationship between Sgas and Ssfr among spiral galax- 
ies, mainly in the Hi-dominated parts. Figure [5] shows 
the relationship between Ssfr and Shi ~ Sgas for 6 H I- 
dominated galaxies. Color contours for the two largest 
galaxies are coded as in Figure [H For these two galaxies 
we plot points from the radial profiles (Figures [2] and [3]) 
on the same plot as black crosses. Because the other 4 
dwarf galaxies are small, we show scatter plots instead 
of density contours. Figure [5] shows individual plots for 
Ho II, IC 2574, NGC 2976 and NGC 4214. Note that the 
remaining galaxies. Ho I, DDO 154, DDO 53, M81 DwA 
and M81 DwB are so small that they yield only 1-10 
sampling points each. We include these data only later 
in Figure [T^l which shows aggregate data for all of our 
dwarf irregular galaxies. 

Figure [5] shows that these galaxies display the same 
saturation value for Shi as the large, centrally H2- 
dominated spirals. This is somewhat surprising, as one 
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Fig. 5. — Contour and scatter plots respectively of EgpH. versus Shi for 2 H l-dominated spirals and 4 dwarf irregular galaxies at 
750 pc resolution. Contour levels, color coding and the dotted diagonal lines of constant SFE are identical to Figure |4] For NGC 925 and 
NGC 2403, crosses show average measurements over tilted rings from the radial profiles. For the 4 dwarf irregulars we show scatter plots 
instead of contours due to the lower number of available sampling points. The dashed vertical line indicates the Sjji saturation level that 
was found for the spirals. 
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would expect conditions in the ISM of many of these 
galaxies to be less favorable to the formation of H2 from 
H I because of comparatively low metallicities (and thus 
lower dust content), shallow potential wells, lower gas 
densities, and more intense radiation fields. One might 
therefore have expected large reservoirs of H i to survive 
in these galaxies at columns where the ISM is mostly 
molecular in a spiral. 

For the most part these galaxies show low SFEs and 
a steep distribution of Ssfr as a function of Shi. This 
includes IC 2574, Holmberg II, and all of the 'small' irreg- 
ulars not plotted. The notable exception is NGC 2976, 
which shows an SFE similar to that found in spiral galax- 
ies and a clear relationship between Ssfr and Sgas- As 
with NGC 4736, this is driven largely by the presence 
of low gas surface densities in NGC 2976 that are rarely 
found in the optical disks of many of the other dwarf 
galaxies. NGC 2976 shows a steadily declining H i profile, 
perhaps curtailed by interactions with other members of 
the M81 group. 

We know from mapping and single dish measurements 
that these galaxies are not H2-dominated (for a Galactic 
conversion factor), so Ehi is likely to be a good proxy 
for the total gas over most lines of sight. However, if the 
CO-to-H2 conversion factor were to vary dramatically, 
as has been suggest e d for dwarf i rregular galaxies ( e.g. , 
iMadden et all 119971 : llsraeil Il997t iLerov et al.l [2007I ). or 
Eh2 > Ehi locally, then Ehi may severely underestimate 
Egas towards the star forming peaks. The topic is too 
complex to address here, but we note the sense of the 
uncertainty: if Egas > Ehi, then the points in Figure 
[5] will move to the right, to higher gas surface densities 
for the same Esfr- Most dwarf irregulars already show 
lower SFEs than most spiral galaxies; the inclusion of 
a substantial reservoir of molecular gas would lower the 
SFE further, thus widening the difference. 

3.7. Dependence on Resolution 

We have so far considered the relationship between gas 
and star formation at 750 pc spatial resolution. Cloud 
formation, stellar feedback, and indeed a breakdown in 
our SFR tracers may all be strong functions of spatial 
scale. Therefore we investigate how the relationships be- 
tween Esfr and Egas or Eh2 change with spatial scale. 
To do this, we convolve all of our data to a variety of 
spatial resolutions (see H2.3P spanning the range from 
the original (highest) resolution to 10 kpc and repeat 
the above analysis. 

Figure [H] shows the distribution of Egas versus EgpR 
for a range of spatial resolutions in two spiral galaxies, 
NGC 5194 (M51) and NGC 6946. We show EgpR versus 
Egas for 4 spatial resolutions, the original (best) spatial 
resolution, and then 500 pc, 750 pc and 1 kpc. Contours 
and other details are identical to Figure [H Aside from 
the expected decrease in the number of independent sam- 
pling points and some narrowing of the distribution as a 
result of averaging, we do not observe any strong effect 
due to spatial resolution. 

To look more quantitatively at the effects of resolu- 
tion on our results, we fit power laws to the data at each 
resolution. As noted before, a single power law is an in- 
adequate description of the data for many galaxies; the 
fits here, however, serve as a shorthand to characterize 
the data distributions. Figure [7] shows obtained power 



law parameters as a function of resolution for 6 of our 7 
spiral galaxies (star formation in NGC 4736 is too con- 
centrated for this exercise) . The fits become more uncer- 
tain as the number of sampling points dwindle and aver- 
aging substantially lowers our dynamic range for many 
galaxies. Therefore we can carry out this test between 
the best available resolution and a spatial resolution of 
1 kpc. The left panels in Figure [7] show the power law 
indices, N, as a function of resolution; the right panels 
show the power law coefficients, A. For the upper panels, 
we fit all CO data above the respective sensitivity limit 
(see § 12.2.211 . For the bottom panels, we fit all total gas 
data within 0.4 r25. 

Figure [7] shows that the derived fits change slowly as 
the spatial resolution decreases from 200 — 500 pc to ^ 
1 kpc. If stellar feedback and cloud formation exert a 
strong infiuence on the relation between star formation 
and neutral gas (and indeed one would expect them to) 
then they do so on scales 5, 300 pc. Th is finding agrees 
with results from iTamburro et al.l ()2008[ ). who measured 
the offset between H i and star forming peaks near spiral 
density waves in some of the same datasets we study. 
They found a typical offset of ~ 100 pc or less between 
regions of star formation and peaks in the H i maps. 

4. COMBINED DISTRIBUTIONS 

In the previous section we saw that the relationship 
between Egas and Esfr varied among spiral galaxies and 
between spirals and Hi-dominated galaxies. We saw 
that the value at which the atomic gas saturates appears 
constant across our entire sample and we found a well- 
defined power law index relating Eh2 and Esfr- In this 
section we collapse the individual contour plots in Figure 
m into combined pixel-by-pixel plots for all of the galax- 
ies in our sample in order to be able to draw general 
conclusions about the star formation law in our galaxy 
sample. 

Figure [8] shows the distribution of Esfr vs. Ehi (top 
left), Eh2 (top right), and Egas (middle right) for all 
sampling points in our 7 spirals. We also assess the im- 
pact of our specific choice of star formation tracer. The 
bottom two panels in Figure [5] show the relationship be- 
tween Esfr and Egas using Ha (bottom left) and, alter- 
nativel y, a combinati o n of H a and 24 /zm emission fol- 
lowing Ic2iittreE^ ()2007D but here applied pixel- by- 
pixel (bottom right). The bottom two plots include data 
for 6 spiral galaxies: NGC 628, NGC 3184, NGC 3521, 
NGC 5055, NGC 5194, and NGC 6946. AU plot pa- 
rameters are identical to those in Figured) In order to 
turn the Ha and the combination of Ha and 24 /im emis- 
sion into SFRs, we assume the same IMF we use for our 
FUV and 24 /im based SFR maps. The Ha-only plot 
furthermore includes a correction for 1.1 magnitudes of 
extinction for every sa mpling point (a typ ical value for 
integrated galaxy disks. lKennicuttlll998b( l. The dotted 
horizontal line in all panels indicates the estimated sensi- 
tivity limits for Esfr - For ease of comparison, the orange 
contour from the middle right plot is overplotted as the 
black contour on both bottom plots. 

One finds that the bottom two distributions in Figure 
[HI which are based on different SF tracers, agree very well 
with the distribution that represents the SF tracer used 
throughout this paper (black contour in bottom panels 
and middle right panel). That is, our choice of star for- 
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Fig. 6. — Contour plots of Ssfr versus Sgas for NGC 5194 (M51) and NGC 6946 at 4 different spatial resolutions: the individual 
maximum resolution, and then 500, 700 and 1000 pc. Contour levels and other features of the plots, including the diagonal dotted lines of 
constant SFE, are identical to those in Figure |4] Although the number of independent data points dwindle, we do not see the distribution 
change markedly with resolution. 




Fig. [6l — continued. 
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Fig. 7. — Power law fit parameters, power law index A'^ and coefficient A (see Equation [2]l , for 6 of our 7 spiral galaxies (NGC 4736 is 
omitted). Left panels: N vs. spatial resolution. Right panels: A vs spatial resolution. The top panels show fits to HERACLES/BIMA 
SONG CO data, the bottom panels to the total gas data within 0.4 r25. Both power law parameters vary only weakly with changing spatial 
resolution. 
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mation tracer does not appear to significantly affect the 
derived relationship between Ssfr and Sgas- 

4.1. Hi Saturation in the Combined Distribution 

The top left panel of Figure [8] shows Ssfr vs. Shi- 
This plot clearly demonstrates the saturation effect dis- 
cussed for the individual galaxies (see ^ 13. 5p and shows 
that it is a universal feature in our sample of spirals. 
The middle left plot shows normalized histograms of 
logEni and logSH2- The Hi shows a clear trunca- 
tion near the threshold. By contrast, the H2 shows 
no such cutoff. This indicates that the H i saturation 
corresponds to a phase transition from an atomic to a 
molecular ISM. The Hi saturation value we quote of 
S_f//, saturation ~ 9 M0 pc^^ represents the 95th percentile 
of the H I distribution. This value coincides w ith th e H i 

l|200l US- 



saturation value derived from IWong fc Blitz 
ing radial profiles in a sample of 6 molecule-rich spiral 
galaxies and agrees with our own radial profiles. 

4.2. Hi, H2, Total Gas, and the Star Formation Law 

The top two panels in Figure [8] show that Ehi and 
Eh2 relate very differently to Ssfr- In the top left plot, a 
narrow range of Sjji corresponds to a large range of Esfr 
and/or SFE. Over less than one order of magnitude in 
Ehi , EsFR covers ~ 3 orders of magnitude and the SFE 
spans 2 orders of magnitude. Ehi therefore cannot be 
used to predict either Esfr or the SFE in spiral galaxies. 
By contrast, the top right panel shows that Eh2 exhibits 
a clear, monotonic relationship with Esfr with a slope 
of approximately unity (see 5 14. 3p down to the sensitivity 
limit of our CO data. 

It is clear from Figure [8] that a single power law can 
describe the relationship between Esfr and Eh2 while 
the distributions of Esfr, and Egas show a clear 'knee' at 
the transition from an H i to an H2-dominated ISM. 

4.3. The Combined Molecular Schmidt Law 

We saw above that the common features in our data 
are the saturation of the H i and the linear relationship 
between Esfr and Eh2- Here we examine this 'molec- 
ular Schmidt law' for the data from all spiral galaxies 
combined. In the top left panel of Figure [9] we show 
the results of a power law fit to Esfr vs. Eh2 for CO 
data from HERACLES/BIMA SONG in spiral galaxies 
(where Egas ~ Eh2)- Contour levels and other plot pa- 
rameters are identical to the ones in Figure |4l The power 
law fit is shown as a solid black line. The best fit param- 
eters are N ^ 1.01 and A = -2.12; the RMS scatter of 
Esfr, about the fit is 0.19 dex, a factor of ^ 1.5. In this 
regime, our data suggest a direct proportionality between 
Esfr, and Eh2 ■ The derived mean H2 gas depletion time 
(including helium) is 2.0 ■ 10^ yrs with an RMS scatter 
of 0.8 • 10^ yrs. 

In the middle left panel the CO data comes from BIMA 
SONG alone. This allows us to plot the data at a slightly 
higher resolution of 500 pc and check the robustness of 
our results to changing CO maps and resolution. All 
other data and plot parameters are identical to the top 
left panel. We obtain an identical slope compared to that 
in the top left panel, N — 1.03 versus N — 1.01, with 
a similar coefficient and scatter. 

The top and middle right panels show Monte Carlo re- 
alizations, where we compute Esfr using Egas from the 



data and the fit parameters quoted in the left panels in- 
cluding the measured lognormal scatter. The agreement 
between the left and right panels provides a qualitative 
check that a single power law is a good description of our 
data in this regime. 

The bottom left and right panels show how the choice 
of star formation tracers affects the derived molecular 
Schmidt law. We use the HERACLES/BIMA SONG 
maps (at 750 pc resolution) and replace our Esfr maps 
with maps derived from Ha emission (bottom left plot) 
and a com bination of Ha and 24 /im emission (bottom 
right plot. ICalzetti et al.ll2007h . The Ha includes a 1.1 
magnitude correction for internal extinction. For these 
data we obtain N = 1.11 and N = 1.18 respectively. 
As in Figure [H the choice of a specific SFR tracer has 
only marginal impact on our results. 

We thus derive a best fit molecular Schmidt law of 

^SFR — tU Zjjj2 . (O) 

for the ensemble, identical to within the uncertainties to 
TV = 1.0 ± 0.1 and A = -2.1 ± 0.2 found for individual 
galaxies (see ij |3.1l and ii l3.2p . The uncertainties of 0.2 
take into account variations in SF tracers, substitution 
of CO maps and scatter in the data. They do not refiect 
variations in the CO-to-H2 conversion factor, the IMF or 
systematics in our methodology. 

5. THE STAR FORMATION LAW AND ENVIRONMENT 

We saw above that there are strong variations in the 
SFE within many spiral galaxies and between spirals and 
H i-dominatcd galaxies. In this section we attempt to 
link variations in the EsFR-Egas relation outside the H2- 
dominated regime in the spirals to variations in environ- 
ment. We approach this in two ways: by measuring how 
the SFE changes with radius within spiral galaxies and 
by comparing Hi-dominated galaxies to spirals. 

5.1. The Radial Dependence of the SFE 

Our data show that below Egas ~ 9 M0 pc~^, the ef- 
ficiency with which gas forms stars varies dramatically 
at a given gas surface density. An obvious explanation 
is that star format i on th r esholds of the sort d is cusse d 
by. e.g.. iKennicuttI (I1989D. iMartin fc KennicuttI (|2001h . 
ISchavel (|2004( l or iLerov et all (|2008af ) are affecting the 
relationship between gas and star formation. We sam- 
ple out to r25 and Martin fc Kennicutt (2001) find that 
the last H 11 region often falls within or near this radius. 
The thresholds described by various authors stem from a 
variety of physics. Shear, Coriolis forces, pressure, metal- 
licity, and passage through spiral arms may all play key 
roles in the formation of gravitationally bound molecular 
clouds, the necessary prerequisite to star formation. 

An exhaustive investigation of all of these other quan- 
tities is beyond the s cope of this paper b ut is the topic of 
a companion paper (jLerov et al.ll2008£l l. However, most 
of the quantities that have been proposed as critical to 
star formation vary radially. Therefore, in Figure [TOl we 
show a scatter plot Esfr vs. Egas for all 7 spirals. Dif- 
ferent colors indicate the galactocentric radius of each 
sampling point. The plotted data are otherwise identical 
to the middle right plot in Figure [H The deprojected 
galactocentric radius of each point, normalized to the 
optical radius of the galaxy, determines the color of the 
point: points within 0.25 r25 are black; sampling points 
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Fig. 8. — Sampling data for all 7 spiral galaxies plotted together. Top left: SsFR vs. Shi; top right: Sgpj^ vs. Sh2; middle right: Sgpj^ 
vs. Sgas- The bottom left and right panels show SgpR vs. Sgas using Ha and a combination of Ho and 24 /zm emission as SF tracers, 
respectively (for a subsample of 6 spirals). The sensitivity limit of each SF tracer is indicated by a horizontal dotted line. The black contour 
in the bottom panels corresponds to the orange contour in the middle right panel and is shown for comparison. The vertical dashed lines 
indicate the value at which Sjji saturates and the vertical dotted lines (top right and middle left panel) represent the sensitivity limit of the 
CO data. The diagonal dotted lines and all other plot parameters are the same as in Figure 3] The middle left panel shows histograms of 
the distributions of H I and H2 surface densities (normalized to the total number of sampling points above the respective sensitivity limit) 
in the sample. 



The Star Formation Law in Nearby Galaxies on Sub-Kpc Scales 



19 







O 




Fig. 9. — The molecular Schmidt law (EgpR versus Sh2) in spiral galaxies. Green, orange and red cells show contours of 1, 2 and 5 data 
points per cell. The black solid line shows the best fit power law for every panel. The power law fit parameters are given in every panel. 
All other lines are the same as in Figures |4] and |8] The top and middle right hand panels show Monte Carlo realizations of the best fit 
power laws using data values for Sgas, observed scatter and best fit parameters. All other panels show data distributions. 
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between 0.25-0.5 are red, those between 0.5-0.75 
are orange; and sampling points from 0.75-1.0 are 
green. 

Figure fTOl shows that in the centers of spirals, where the 
gas is mostly molecular, the H2 depletion time is nearly 
constant (see S I4.3p . At larger radii, where the ISM is 
dominated by Hi, the SEE of gas at a particular sur- 
face density depends on the galactocentric radius: The 
further out in the disk, the lower the SEE. 

Eigure [11] separates the data from Eigure [TO] into in- 
dividual plots of SsFR versus Sgas for specific ranges of 
galactocentric radius. We use the same radius bins as in 
Eigure [TO] and show a separate contour plot for each bin. 
The upper left panel shows the black sampling points 
(< 0.25 r25) from Figure [TO] the upper right plot shows 
the red points (0.25-0.5 r25), the lower left plot the or- 
ange (0.5-0.75 7-25), and the lower right plot shows the 
green sampling points (0.75-1.0 r25). In all plots, we 
quote the mean SEE in that radius bin and its RMS 
scatter. 

Eigure [TT] shows that where Shi > 5]h2, ^sfr is a 
function not only of Sgas, but also of local conditions 
that vary strongly with galactocentric radius. These in- 
clude, e.g., metallicity, stellar surfa ce density, gas pres- 
sure, galactic rotation, and shear. iLerov et al.l (|2008ar ) 
explore the relationship between these quantities, Ssfrj 
and Sgas using the same data we use here. 

5.2. Hi-dominated Galaxies and the Outer Disks of 
Spirals 

We saw in ii l3.6l that the SEE in H i-dominated galaxies 
is lower than that in spiral galaxies. In the previous sec- 
tion, we showed that the SEE varies dramatically within 
a spiral galaxy as a function of radius. Here we compare 
these two findings. 

Eigure [TO] shows the combined distribution of Ssfr 
versus Shi for all dwarf irregular galaxies in our sample 
(i.e., the Hi-dominated galaxies from Table [T] omitting 
NGC 925 and NGC 2403 that would otherwise dominate 
the distribution). Wc show the same data for the Hi- 
dominated galaxies in all 4 panels. In each panel, we 
overplot a black contour that shows Ssfr versus Sgas 
for the spiral galaxies from a particular radial bin. This 
black contour in each panel of Eigure [TO] corresponds to 
the lowest (green) contour of the respective panel in Eig- 
ure [H 

One finds that the distribution of Ssfr versus Sgas for 
the Hi-dominated/dwarf galaxies overlaps the distribu- 
tion of the outer disks of the spiral galaxies (compare 
bottom right panel of Eigure [TO]) . Many conditions are 
similar in these two regimes - low metallicities, low dust 
to gas ratios, high atomic gas fractions, and compara- 
tively weak stellar potential wells. This plot suggests 
that these shared environmental factors lead to a simi- 
lar relationship between gas and star formation in both 
regimes. 

We note that one environmental fact that dwarf irreg- 
ulars do not share with the outer disks of spiral galax- 
ies is differential rotation. Dwarf galaxies tend to have 
nearly solid body rotation curves and correspondingly 
low shear, whereas rotation curves tend to be flat in the 
outer parts of spirals. This implies that shear alone may 
not be the driving force regulating the SEE. 



5.3. The Molecular-to- Atomic Gas Ratio Sh2 /Sri as a 
Function of Radius 

In ^ I5.1l we saw that where Shi > Sh2, the SEE varies 
strongly with radius. In combination with our finding 
that Ssfr oc Sh2 (see ij |4.3p . this leads us to expect that 
the H2-to-Hi ratio also varies strongly with radius. The 
right panel in Eigure [13] shows Sh2/Shi as a function of 
normalized galactocentric radius for all 7 spiral galaxies. 

We find that the H2-to-H i ratio decreases as a function 
of radius. But the sensitivity of our CO data limits the 
measurement beyond ~0.5r25. An alternate approach is 
therefore to extrapolate Sh2 from Ssfr using the direct 
proportionality between the two that we have established 
in S I4.3I for the inner part of the optical disk. The EUV 
and 24 /im data are more sensitive than our CO maps 
and so offer more significant measurements in the outer 
disk. Eor this approach, we assume that the relationship 
that we fit between Ssfr and Sh2 holds across the disks 
of our spiral galaxies and infer the ratio Sh2 (Ssfr) /^hi 
for each sampling point. This is shown in the left panel. 
Eor comparison, the green contour in the left panel is 
overplotted in the right panel as a black contour. 

One finds that the two approaches agree reasonably 
well; the ratio measured from CO overlaps the ratios 
inferred from Ssfr- Both panels show that the phase 
of the ISM is a well-defined function of local conditions 
that vary with radius. A detailed investigation is be- 
yond the scope of this paper but we explore this topic 
(o ften referred to as 's tar formation thresholds') further 
in iLerov et al.l ((20083) ■ Here we show that an exponen- 
tial fit, which is shown by the dashed black line in both 
panels, can provide a good description of the relationship 
between H2 , H i and the SER where Sri > Sh2 • The 
exponential fit yields a scale length of 0.2 r25 in both pan- 
els. The radius where Sh2 = 5^hi can be derived from 
the coefficients of the fits. The result is identical in both 
cases: 0.43 7*25 in the left panel and 0.46 r25 in the right 
panel respectively. 

6. SUMMARY & DISCUSSION 
6.1. Summary of Our Work 

Using data from THINGS, HERACLES, BIMA SONG, 
the GALEX NGS and SINGS, we derive Sri, Sh2, and 
Ssfr at 750 pc resolution across the optical disks of 7 
nearby spiral galaxies with H2-dominated centers and 11 
H i-dominated late-type galaxies. We use these datasets 
to make the first pixel-by-pixel analysis of the star for- 
mation law in a significant sample of nearby galaxies. 

We find two relationships common throughout our 
data. Eirst, a molecular Schmidt law with index N = 
1.0 ± 0.2 relates Sh2 to Ssfr in our sample of spirals. 
This may also be described as a total gas Schmidt law in- 
side ^ 0.5 r 25, where the ISM of all of our spiral galaxies 
is H2-dominated. The average molecular gas (including 
helium) depletion time is 2.0 • 10^ years with an RMS 
scatter of 0.8 • 10^ years. This relationship holds for indi- 
vidual galaxies as well as for the combined distribution 
and is robust to substituting BIMA SONG for HERA- 
CLES data, using different SER tracers, or changing the 
resolution from ~ 300 pc to ^ 1 kpc. 

The second common feature of our data is that Shi sat- 
urates at a surface density of 9 M0 pc~^; gas in excess 
of this value is found in the molecular phase in the spirals. 
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Fig. 10.— SsFR versus Sgas for the spiral galaxies with data colored according to galactocentric radius. Every data point represents 
one sampling point; the data are the same as in the middle right plot of Figure |8] We show the normalized galactocentric radius of every 
sampling point via its color, where, in units of r25, the colors correspond to: black < 0.25; red 0.25—0.5; orange 0.5-0.75; and green > 0.75. 
The diagonal dotted lines and all other plot parameters are the same as in Figure |4] The data clearly break up according to galactocentric 
radius, with low SFE/high gas depletion time points corresponding to the outer parts of the spirals. 



This saturation is common to spiral and Hi-dominated 
galaxies. This is somewhat surprising because conditions 
in H i-dominated galaxies (such as dwarf galaxies) should 
be less favorable to the formation of H2, which may lead 
one to expect a higher saturation value; a situation that 
is not observed. 

We do not observe a universal relationship between 
total gas surface density and Ssfr- Outside the H2- 
dominated region, i.e., at r > 0.5 r25, the relationship 
between gas and star formation varies both within and 
among galaxies. In some cases, a single power law relates 
total gas and SFR over many orders of magnitudes in gas 
surface density. In other cases, we find a wide range of 
star formation rates at almost the same gas column. As 
a result of this variation, our best-fit power law index, 
N, for the total gas in spiral galaxies ranges from 1.1 to 
2.7. This agrees well with the range of indices found in 
the literature, but does not hint at a universal total gas 



Schmidt law. 

We describe variations in the relation between the 
SFR (Ssfr) and the total gas (Sgas) in terms of the 
star formation efficiency (SFE), i.e., star formation rate 
per unit gas mass (EsFR/Sgas), and the Il2-to-Hi ratio 
(I]h2/Shi)- We show that outside of the central, H2- 
dominated regime of spiral galaxies, the SFE has a strong 
gradient with radius, where the highest efficiency points 
come from the inner disk, and the lowest efficiency points 
are at larger radii. We also show that the pixel-by-pixel 
distributions of late-type, H i-dominated galaxies over- 
lap those of the outer disks of the spirals in our sample 
in EgFR-Sgas parameter space. This implies that simi- 
lar conditions, i.e., low metallicities, weak potential wells 
and low dust content, might drive the SFE in this regime. 

We argue that these observations show a clear link be- 
tween environment and the relationship between gas and 
star formation. We suggest the following scenario: the 
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Fig. 11. — The variation of SgpR versus Sgas with radius in spiral galaxies. Levels of 1, 2, 5 and 10 points per cell are shown as 
green, orange, red and magenta contours. The diagonal dotted lines and all other plot parameters are the same as in Figure |4] Every 
plot represents only sampling data from a certain range in normalized galactocentric radius, corresponding to a particular color of points 
in Figure [TOl Upper left panel: < 0.25 r25 (black points); upper right: 0.25—0.5 r25 (red points); lower left: 0.5—0.75 r25 (orange points); 
lower right: 0.75-1.0 r25 (green points). Data from outer galaxy disks display lower SgpR for the same Sgas compared to data from the 
inner disks. 



observational star formation 'law', as it is now under- 
stood (e.g., K98), is a molecular phenomenon. The tran- 
sition from Hi to H2 and the subsequent formation of 
stars is not purely a function of the total gas surface 
density. Instead, other physics sets the ratio of Hi to 
H2. The critical quantity for these processes is a strong 
function of radius and appears to be common to both 
dwarf irregular galaxies and the outer regions of spiral 
disks. This agrees quantitatively with the findings of 
IWong fc Blitd ()2002f ) and qualitatively with t he ideas of 
star formation thresholds d i scusse d by, e.g., ^ Kcnnicutll 
(fT989l) . [Martin fc KennicuttI (|200lh . iLerov et al . (2008^) 
and many others. 

We test these results in several ways and find that 
our conclusions are robust against variations in the star 
formation rate tracer used or the applied spatial reso- 
lution. Moreover, we find the same results when per- 
forming a pixel-by-pixel analysis and when working with 



azimuthally averaged radial profiles. 

6.2. Comparison with Previous Measurements 

As we discussed in §[1] many different studies have ex- 
amined the SF law in nearby galaxies and derived a vari- 
ety of power law indices. These studies have used a wide 
range of SFR tracers, methodologies (disk averaged mea- 
surements, radial profiles or individual apertures) and 
assumptions (e.g., adopted IMF or CO-to-H2 conversion 
factor) . Here we compile the results from several studies 
whose galaxy samples overlap our own to compare the 
data distributions in EsFR-Sgas phase space once differ- 
ent methodologies are accounted for. 

In Figure [HI we therefore compare our Ssfr vs. Egas 
pixel-by-pixel distributions and our radial profile data 
to measurements from other authors after matching the 
assumptions which have gone into each of their analyses 
to ours. In each panel, the colored contours represent the 
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Fig. 12. — A comparison between H l-dominatod/dwarf irregular galaxies and different radial regimes of spiral galaxies. All four panels 
show SsFR versus Ejji for the dwarfs in colored contours. Green, orange, and red contours show 1, 2, and 5 sampling points per cell. 
The diagonal dotted lines and all other plot parameters are the same as in Figure |4] Overplotted is the lowest (green) contour from the 
corresponding panel in Figure [TT] Thus each panel compares the distribution of data from the dwarfs to that in the spiral galaxies from a 
particular radial range. The best agreement is seen in the bottom right panel, in which the black contour shows data from 0.75-1.0 r25 in 
spiral galaxies. 



pixel-by-pixel data from our study; these distributions 
are identical to the ones in Figure [H The black crosses 
represent the datapoints from our radial profiles (Figure 
[5]) and are also identical to the crosses in Figure I?) 

The black dots in Panel A and the black circles in 
the other panels show data from previous papers. These 
data were adjusted to match our assumptions regarding 
the adopted IMF (see ^ I2.3.2p . the inclination, the CO-to- 
H2 conversion factor, and the lco{J = 2 ^ ^)l^co{J = 
1 — > 0) line ratio. We mention the adjustments that were 
made for each dataset in the discussion of the individual 
panels of F igure [Til below. 

Panel A: iKennicutt et al.l ()2007l K07) derive a power 
law index of 1.37±0.03 relating molecular gas and SFR 
surface densities in M51 (and an index 1.56±0.04 for 
the total gas). They measure Egas and Ssfr by plac- 
ing 520 pc apertures on Ha and 24 /.tm emission peaks 



(black dots in Panel A in Figure [H]) . K07 use Paa and a 
combination of 24 /im and Ha emission to estimate Ssfr- 
They use the same Hi and CO data we do. 

Shown (Figure [U A) is a comparison for NGC 5194 
(M51) between our data and the measurements from 
K07, which we have adjusted to match our CO-to-H2 
conversion factor and IMF. One finds that our data 
distribution generally agrees quite well with that from 
K07. Nevertheless, the regime that both studies probe is 
slightly different: whereas we sample the entire optical 
disk, K07 focus on apertures mainly in the spiral arms 
of M51. As a consequence, our data distribution extends 
to lower Egas while their distribution emphasizes slightly 
higher Esfr and Egas- 

The power law slope of 1.37±0.03 derived in K07 for 
the molecular gas is steeper than the value of 0.84 that we 
derive for M51 (see Figure |4] and Table [2|). The steeper 
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Fig. 13. — The ratio of molecular to atomic gas, Sh2 / Shi, a function of galactocentric radius for the spiral galaxies. In the left 
panel, Sh2 is inferred from EgpR using the direct proportionality between the two that we have established in ii l4.3l In the right panel, Sh2 
from HERACLES is used directly. The black contour in the right panel is identical to the green contour in the left panel and is shown for 
comparison. Green, orange, and red cells show contours of 2, 5 and 10 sampling points per cell (cell sizes: Ax = 0.02 r25. Ay = 0.08 dex). 
The filled black circles show the median Ssfr per 0.1 r25 radial bin, the error bars represent the la scatter in each bin. We carry out 
exponential fits to both distributions (dashed black lines). 



slope in K07 is driven mainly by the fact that they mea- 
sure low Ssfr for some apertures at log(I]gas) between 
~1 and 1.8 (compare black points in panel A). This dif- 
ference arises because K07 subtract a local Ha and 24 /im 
background for each aperture. This has marginal impact 
at high Ssfr, but significantly affects low Ssfr measure- 
ments where the contrast with the local background is 
low. This comparison shows that despite using the same 
gas data and similar SFR tracers, the applied methodol- 
ogy for deriving SFRs and the specific choice of sampling 
method have non-negligible impact on the derived power 
law slopes. 

In Panel B we compare ou r M51 distribution to data 
adopted froml Schuster et"an ()2007[ ). which are adjusted 
to match our IcoiJ = 2 — > 1) /lco{J =1^0) line ratio, 
IMF and C0-to-H2 conversion factor (we note that they 
apply a conver sion factor equ a l to o ne quarter of our 
adopted value). ISchuster et al.l ()2007D use radial profiles 
and obtain a power law index of 1.4±0.6 when fitting 
SFR and total gas surface densities in M51. They use 
radio continuum emission at 20 cm to derive Ssfr- 

In this panel, we plot their data after adjusting for 
differences in the applied line ratio, IMF and CO-to-Il2 
conversion factor and find good agreement for the distri- 
butions. The dis crepancy at sma l l Sgas arises because the 
Hi data used in ISchuster et al.l (f2007) are less sensitive 
than our map. This leads to decreasing Hi columns at 
larger radii in their radial profile while our profile remains 
flat. T he central radial profile point from ISchuster et al.l 
()2007[) shows a particularly high Ssfr- This may be due 
to the AGN in M51, which is a strong source of radio 
continuum emission. 

Panel C and D: Both panels show radial profile 
data adopted from IWong fc BliTg (|2002D . panel C for 
NGC4736 and pa nel D for NGC5055 . In both pan- 
els, the data from iWong fc Blit3 (|2002l ) are adjusted to 
match our IMF and adopted inclinations. They derive 
their SFRs from Ha emission and for their plotted data- 



points in Panels C and D, we have adopted their constant 
global extinction c orrection of 1.5 mag. We note that the 
radial profiles from I Wong fc Bli"t3 (|2002f ) only extend out 
to ~150" and -300" for NGC4736 and NGC5055, re- 
spectively, while ours reach the edge of the optical disks 
(~230" and ~350", respectively). 

For NGC4 736, the observed Ssfr for the 
IWong fc Blitl (|200l profile are lower than those mea- 
sured by us (apart from their 'bump' at log(Sgas) « 1-5). 
For this galaxy, their derived Ssfr depends sensitively 
on the method used to correct the Ha emission for 
extinction: mid-IR (24 /im), fixed extinction, and a 
gas column based extinction estimate yield results that 
differ by a factor of about two. In the case of NGC 5055 
(Panel D), the agreement between both distributions is 
excellent for all Sgas- 

Pa nel E compares dat a for NGC 6946 adopted 
from ICrosthwaite fc Turneil (|2007f ). adjusted to match 
our adopted inclina t ion and IMF, to our data. 
ICrosthwaite fc Turned ()2007f) use radial profiles and find 
a proportionality between the molecular gas and the SFR 
surface density (implying a power law slope of ~1 for the 
molecular gas) . They derive their SFRs from FIR and ra- 
dio continuum emission. We will, however, only use the 
FIR based SFRs in our comparison, because it is easier 
to match assumptions for this case. Our radial profile 
extends to larger radii (and thus includes lower Sgas val- 
ues) and because ICrosthwaite fc Turner] ()2007( ) average 
emission in 30" wide rings (our study: 10"), they have 
fewer datapoints and their Ssfr and Sgas show less dy- 
namical range. Despite an offset in Ssfr, likely due to 
the different SFR tracers, both distributions are in good 
qualitative agreement. 

Panel F: We compare pixel-by-pixel data from 
NGC 3351, which is n ot part of our curr ent sample, to 
inner disk data from iJogee et al.l (|2005f ). Their data 
are adjusted to match our IMF, CO-to-H2 factor and 
adopted inclination (41°). They use radial profiles to 
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Fig. 14. — Sgpfi versus Egas for individual galaxies from this study and from previous analyses. The colored contours represent the 
pixel-by-pixel data for spiral galaxies from this study; the distributions are identical to the respective ones in Figure|4] The diagonal dotted 
lines and all other plot parameters are the same as in Figure|4] Overplotted are datapoints from our radial profiles (compare Figures [2] and 
|4]l as black crosses. We compare our data to measurements from other analyses, which were carried out either by using individual apertures 
(panel A; each black dot represents the measurement within an individual aperture) or radial profiles (other panels; each filled black circle 
represents a radial profile data point). The references are given in each panel. We adjust these data to match our adopted inclinations, 
IMF, CO line ratio and CO-to-H2 conversion factor. In general, there is good agreement between our data and the measurements from 
other studies, but the comparison also shows that the distributions are quite sensitive to some of the underlying assumptions, e.g., how 
one corrects for extinction or which CO-to-H2 conversion factor is applied. 
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sample the inner Ikpc of NGC3351. Because they re- 
strict their analysis to the inner disk, they only consider 
molecular gas. They derive their SFRs from Br7 emis- 
sion. This comparison allows us to study a regime of 
relatively high Ssfr and Sgas- Due to the coarser res- 
olution of o u r dat a (750 pc) compared to the data from 
iJogee et al.l ()2005l . 120pc for NGC3351), we can only 
compare 3 radial profile points (black crosses) and a few 
individual pixel-by-pixel sampling poi nts (highlighted as 
black circles) to the measurements of iJogee et al.l ()2005l 
black filled circles). We find an offset in terms of Ssfr, 
probably due to the different SFR tracers used. Both 
distributions show good qualitative agreement neverthe- 
less. 

These comparisons show that once the different as- 
sumptions of all measurements are accounted for, the 
data distributions from the different studies agree well 
in general. This appears to be relatively independent 
of the specific SFR tracer or the samphng method (e.g., 
pixel-by-pixel or radial profiles) that is applied. The res- 
olution of the data seems to matter only marginally, as, 
even when disk averaged measurements are included, the 
data still populate the same regime in SsFR-Sgas phase 
space (as we will discuss in ij |6.3p . 

When one corrects for different assumptions, the CO- 
to-H2 conversion factor is especially important and mat- 
ters particularly in the inner part of galaxies, where the 
molecular gas dominates the total gas budget. The fac- 
tors originally applied in the above studies differ by as 
much as a factor of 6. If one does not correct for different 
conversion factors, the high Sgas end of the distributions 
may change quite drastically, which can significantly in- 
fluence derived power law slopes. 

Nuclear activity, i.e., the presence of AGN, may also 
play a role for radial profile datapoints from the centers 
of galaxies, depending on the applied SFR tracers (see 
discussion for Panel B). This is only a minor issue when 
sampling pixel- by-pixel, where only very few pixels from 
the center contribute. 

6.3. The Molecular Schmidt Law In Various Regimes 

Figure [T5] shows our combined pixel- by-pixel distribu- 
tion together with all the literature measurements dis- 
cussed in SI6.21 m easurements from starburst galaxies 
(|Kennicuttj|1998af) . and data from low surface brightness 
galaxies (LSB, Wyder et al., in prep.). Our data and 
those from the literature sweep out a clear, consistent 
trend in SFR-gas space. The 'kink' in the pixel-by-pixel 
distribution (indicated by the left vertical line) , refiecting 
the saturation of H i at Sgas ~ 9 pc~^ as discussed 
in §[33] and § 14.11 is also visible in the data distribution 
from the other studies. 

In the following we will discuss how the results from 
this study relate to the work of K98, which is the only 
study so far to explore the SF law over the full range of 
SsFR and Egas shown in Figure I15[ and we will discuss 
a possible transition that may be expected at around 
Egas ~ 200 M0 pc^^ (right vertical fine in Figure fTS]) . 

In this paper, we find Esfr oc H]^"^^'^, whereas K98 
found a power law with slope N — 1.40 ± 0.15 relating 
EsFR and Egas- The fit of K98 depends on the contrast 
between normal spirals, Eh2 ~ 20 pc~^, and high 
surface density starbursts, Eh2 ~ 1000 Mq pc~^. A 



power law index N « 1.5 relating SFR to CO emission 
has been well-established in starbursts at low and high 
redsh i fts by a number of a uthors (e.g.. iGao fc SolomonI 
120041: iRiechers et al.|[2007l ). There may be reasons to 
expect different values of N in starburst environments 
and in our data. Starburst galaxies have average sur- 
face densities far in e xcess of a Galactic GMC (e.g., 
IGao fc Soloi^ [200l [Rosolowskv fc Bliti [2005h . We 
have no such regions in our own sample, instead we 
make our measurements in the regime where Eh2 = 
3 — 50 M0 pc"^. In starbursts, the changes in molecular 
surface density must refiect real changes in the physical 
conditions being observed. 

In our data, Eh2 is likely to be a measure of the fill- 
ing factor of GMCs rather than real variations in sur- 
face density. On the one hand, for our resolution (750 
pc) and sensitivity (Eh2 = 3 Mq pc^^) the minimum 
mass we can detect along a line of sight is ~ 1.5 • 10^ 
Mq. Most of the mass in Galactic GMCs is in clouds 
with Mh2 « 5 • 10^ - 10^ Mq (e.g., lBlit3 Il993h . Con- 
sequently, wherever we detect H2 we expect at least a 
few GMCs in our beam. On the other hand, most of our 
data have Eh2 ^ 50 M0 pc~^. The typi cal surface den- 
sity o f a Galactic GMC is 170 Mq pc"^ (jSolomon et all 
|1987() . These surface densities are much lower than those 
observed in starbursts and are consistent with Galactic 
GMCs filhng < 1/3 of the beam. If GMC properties are 
the same in all spirals in our sample, then for this range 
of surface densities we expect a power law index of = 1 
as Eh2 just represents the beam filling fraction of GMCs. 
Averaging over at least a few clouds may wash out cloud- 
cloud variations in the SFE. A test of this interpretation 
is to measure GMC properties in a wide sample of spirals. 
We note that Local Group spirals display similar scaling 
relations and cloud mass distribution functions so that it 
is hard to distinguish G MCs in M 31 or M 33 from those 
in th e Milky Way (e.g.. lBlitz et al.l [20071 : IBolatto et all 
l2008f) . If this holds for all spirals, then we may indeed 
expect = 1 whenever GMCs represent the dominant 
mode of star formation. The next generation of mm- 
arrays should soon be able to measure GMC properties 
beyond the Local Group and shed light on this topic. 
In that sense, our measurement of A^ = 1.0 ± 0.2 repre- 
sents a prediction that GMC properties are more or less 
universal in nearby spiral galaxies. 

For our results to be consistent with those from star- 
bursts, the slope must steepen near Eh2 ~ 200 M© pc^^. 
This might be expected on both observational and phys- 
ical grounds. CO is optically thick at the surfaces of 
molecular clouds. Therefore, as the filling fraction of such 
clouds for a given telescope beam approaches unity, CO 
will become an increasingly poor measure of the true Eh2 
because of the optical thickness of the CO emission. Even 
if such clouds have Galactic SFEs, the observed relation- 
ship between CO and an optically thin SFR tracer, e.g., 
far infrared (FIR) emission, will steepen. It is also likely 
that physical conditions inside the molecular gas change 
as surface densities exceed that of a typical GMC (e.g., 
iPadoan et al.ll2007l ). If the density increases, the free 
fall time within the gas will decrease^ possibly leading to 
more efficient star formation (e .g., [K rumholz fc McKei 
'2005"). We refer the reader to Krumho lz fc ThompsonI 
(2007 ) for a full treatment of theoretical expectations for 
changing Schmidt laws. 
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Fig. 15. — SsFR versus Egas from this paper in colored contours (compare middle right panel of Figure IS} and for individual galaxies 
from other analyses (see Figure [T4t . The diagonal dotted lines and all other plot par ameters are the same as in Figure |4] Ovcrplotted as 
black dots are data from measurements in individual apertures in M51 IjKennicutt et al. 2007). Datapoints from radial profiles from M51 
ISchuster et al. 2007), NGC 4736 and NGC 5055 ( Wong & Blitz 2002) and from NGC 6946 ( Crosthwaite & Turner 2007) arc shown as black 
filled circles. Furthermore, we show disk averaged measurements from 61 normal spiral galaxies (filled grey stars) and 36 starburst galaxies 
(triangles) from K98. The black filled diamonds show global measurements from 20 low surface brightness (LSB) galaxies (Wyder et al., 
in prop.). Data from other authors were adjusted to match our assumptions on the underlying IMF, CO line ratio, CO-to-H2 conversion 
factor and galaxy inclinations where applicable. One finds good qualitative agreement between our data and the measurements from other 
studies despite a variety of applied star formation rate tracers. This combined data distribution is indicative of 3 distinctly different regimes 
(indicated by the vertical lines) for the star formation law (see discussion in the text). 

If this sketch is correct, then care must be taken re- 
lating the SFR to CO measurements. The efficiency 
with which gas forms stars will depend on what regime 
one considers. Perhaps more excitingly, the reverse is 
also true: by measuring the SFE one can place a re- 
gion or galaxy in cither the 'starburst' or 'GMC/disk' 
regime. This offers the intriguing prospect, for exam- 
ple, of diagnosing the dominant mode of star formation 
in surveys of star formation and molecular gas at high 
redshift that will be a major component of ALMA sci- 
ence. Even with the current generation of millimeter- 
wave t elescopes such coni parisons are possible. For ex- 
ample, iDaddi et al.l (|2008f ) recently showed that molecu- 
lar gas in 'BzK'-selected galaxies at z = 1.5 has the same 



In fact, iGao fc Solomonl ([2f)0l may observe the tran- 
sition directly; in their Figures lb and 3 they show that 
FIR emission (i.e., the SFR) scales linearly with CO 
emission below Mh2 ~ 10^" M0, i.e., for their sample of 
'normal' spiral galaxies. Above 10^° Mq, i.e., for LIRGs 
and ULIRGs, they observe a steeper power law index re- 
lating SFR to CO (while the relationship between SFR 
and dense g traced by HCN emission, remains lin- 
ear). This value, Mh2 ~ 10^° Mq, is about the upper 
limit of molecular gas masses observed in our sa mple of 
spiral galaxies (see Table 5 in lLeroy et al.l[2008al ). Based 
on our results, we suggest that the former (linear) regime 
corresponds to star formation organized into normal spi- 
ral GMCs. 
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FIR-to-CO ratio (and thus presumably SFE) as local spi- 
ral galaxies. For comparison, sub-millimeter ('SCUBA') 
galaxies display dramatically higher SFEs, as t raced by 
their FIR-to-CO ratios fe.g. JCreve et al.ll2005[) . A sug- 
gestive link to this work is that the BzK galaxies show 
signs of larger stellar disks and more extended CO emis- 
sion, completely consistent with star formation proceed- 
ing mostly in a disk populated by GMCs analogous to 
those in nearby galaxies. 
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